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In this paper, we consider maximizing the sum-rate in the uplink of a multi-cell OFDMA network. 
The problem has a non-convex combinatorial structure and is known to be NP hard. Due to the inherent 
complexity of implementing the optimal solution, firstly, we derive an upper and lower bound to the 
optimal average network throughput. Moreover, we investigate the performance of a near optimal single 
cell resource allocation scheme in the presence of ICI which leads to another easily computable lower 
bound. We then develop a centralized sub-optimal scheme that is composed of a geometric programming 
based power control phase in conjunction with an iterative subcarrier allocation phase. Although, the 
scheme is computationally complex, it provides an effective benchmark for low complexity schemes 
even without the power control phase. Finally, we propose less complex centralized and distributed 
schemes that are well-suited for practical scenarios. The computational complexity of all schemes is 
analyzed and performance is compared through simulations. Simulation results demonstrate that the 
proposed low complexity schemes can achieve comparable performance to the centralized sub-optimal 
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scheme in various scenarios. Moreover, comparisons with the upper and lower bounds provide insight 
on the performance gap between the proposed schemes and the optimal solution. 

I. Introduction 

Dynamic resource allocation plays a central role in the air interface design of state-of-the- 
art OFDMA-based cellular technologies. In this paper, we focus our attention on maximizing 
the overall network throughput by optimizing the allocation of resources (i.e., subcarriers and 
powers) jointly in a multi-cell uplink OFDMA network. The goal is to develop efficient resource 
allocation schemes that takes into account the inter-cell interference (ICI) while considering 
universal frequency reuse. The solution of such problem is difficult to achieve optimally due to 
its NP hard combinatorial nature and high dimensionality. 

The sum rate maximization problem is extensively studied for the downlink in OFDMA 
networks. The optimal strategy in the downlink is to separately optimize subcarrier and power 
allocation, i.e., allocate a subcarrier to the user with best channel and then perform water-filling 
over the allocated subcarriers [1]. However, the problem becomes more challenging in the uplink 
scenario due to the individual power constraint at each user. Simply allocating a subcarrier to 
the user with best channel quality may affect the network performance considerably, as some 
active users may have better channel gains but low transmission powers on a specific subcarrier. 

Most of the recent work in the context of multi-cell OFDMA networks (7]|, flU, aims at 
minimizing the overall transmitted power, i.e., linear objective with pre-defined rate constraints. 
In BH, the authors investigated scaling laws for upper and lower bounds of the downlink capacity 
in the asymptotic regime. Furthermore, in some recent literature, low complexity distributed 
game theoretic solutions are also studied. However, the schemes are iterative and optimality is 
not guaranteed [fTOTl . An auction based approach is discussed in ifTTI . where the authors proposed 
a joint auction and dual decomposition based technique for the resource allocation problem. The 
technique is asymptotically optimal as the number of subcarriers in every cell goes to infinity. 
However, this may not be true for finite number of carriers. In summary, all these approaches 
are sub-optimal and no criteria are used to calibrate their performance gap with respect to the 
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optimal solution. 

Motivated by the above discussed facts, we consider the problem of optimized resource 
allocation in the uplink of multi-cell OFDM A networks. Firstly, we compute an upper bound 
(UB) and lower bound (LB) to the optimal average network throughput. Also, we study the effect 
of ICI on the performance of the near-optimal single cell resource allocation scheme proposed 
in [O which leads to another simple lower bound. Simulation results show that this lower bound 
is slightly loose but can be computed easily. Since the computation of the optimal solution 
is exhaustive, we then propose a centralized sub-optimal resource allocation scheme which 
uses a geometric programming (GP) based power control phase in conjunction with a heuristic 
subcarrier allocation phase. The proposed scheme possesses an iterative and computationally 
intensive subcarrier allocation phase. However, it can serve as an effective benchmark for the 
less complex schemes even without the power control phase. Furthermore, the power control 
phase is discussed in this paper for both high and general signal to interference plus noise 
ratio (SINR) regimes. Finally, we propose and evaluate less complex centralized and distributed 
schemes that are suitable for practical implementation. 

The rest of the paper is organized as follows: In Section II, the system model is defined and 
problem is formulated. In Section III, the bounds are derived and their complexity is analyzed. 
In Section IV and Section V, the proposed centralized and distributed schemes are explained. 
Section VI presents numerical results followed by concluding remarks in Section VII. 
Notation: Throughout the paper, we denote the sets of real and complex vectors of N elements 
by ~R N and C^, respectively. Matrices are represented using boldface upper case letters while 
bold face lower case letters are used for vectors. A/"(0, a 2 ) represents a zero mean Gaussian 
random variable with variance a 2 . 

II. System Model and Problem Formulation 

A network of L cells with a set of K users in each cell / is considered. Full reuse of the 
spectrum is assumed in all the cells (i.e., frequency reuse =1). Each base station (BS) is assumed 
to have N orthogonal subcarriers, and each subcarrier can be allocated to a single user per cell. 
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The average network throughput C is a function of both subcarrier and power allocation variables. 
The sum rate maximization problem is formulated as follows using the standard Shannon capacity 
formula, C n ,k,i = log2(l + 1n,k,i), where C n> k,i and 7^*^ represent the throughput and SINR of 
the k th user at n th subcarrier in cell /, respectively: 

L K N 



sr^sr^sr^ , {-, , Pn,k,ih n ,k,i\ 
maximize ^ ^ ^ a nAl log 2 f 1 + j 



(1) 



N 



SUbjeCt tO ^P„,fc,Z < Pfe,max, V/c,V/ (2) 



n=l 
K 



^a nA i = 1, Vn,V/ (3) 



fc=i 



a„, M e{0,l}, Vn,V/,VA; (4) 



In ©, 4 )Z = Yl^=i,j^iYlk=i a n,k,jPn,k,j9n,k,ji represents the cumulative interference at - 
subcarrier in cell / from the users in all other cells, p n k,i denotes the power transmitted by k th 
user at the n th subcarrier in cell /, a n> k,i represents the allocation of k th user at the n th subcarrier 
in cell / and h n k i is the channel gain of k th user at the n th subcarrier in cell I. Constraint © 
implies that the power spent by k th user on its allocated subcarriers cannot exceed the maximum 
available power, Pfc max . For each cell, we collect the power allocation variables p n ,k,i in a vector 
Pn,« = [Pn,i,hPn,2,h ■■■■,Pn,K,i] an d then stack all the vectors in a power matrix P/ of cell I where 
Pz G M. NxK . Constraint ([3]) restricts the allocation of each subcarrier to only one user. The 
channel gains h n> k,i and binary allocation variables a n ^,i are stacked up similarly in the matrices 
H/ and A/, respectively, where Aj, H; E ~R NxK . Moreover, we define g n ,k,ij as the interfering 
gain from the k th user in cell I to cell j, Vj 7^ I at n th subcarrier. We collect these interfering 
gains into a vector g n /j = [g n ,i,ij, g n ,2,ij 9n,K,ij] and then stack all the vectors in a matrix 
Gij e R NxK . 



n 
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A. Optimal Problem Formulation in High SINR Regime 

Assuming perfect knowledge of channel gains at a centralized controller, the optimal solution 
for (OQ) can be computed in the high SINR regime by an exhaustive search over all possible 
combinations of the allocations. For each possible allocation, optimum powers can be computed 
by transforming (OQ) into a GP. Note that the power allocation problem is in itself a known 
non-convex problem for the general SINR regime lfl4ll . However, in the high SINR regime the 
problem becomes a convex GP problem. For a given set of allocation variables and considering 
a high SINR regime, the objective function in (Q~|) can be rewritten as follows: 

L K N f h \ 

maximize £ E ^TfT ) (5) 

Pn ' kd 1=1 k=l n=l V ° + J n,Z / 

Maximizing the SINRs is equivalent to minimizing the interference to signal ratio: 

minimize V] Y] Y] a n k jlog 2 ( 7-^ J (6) 

tttttl \Pn,k,lhn,k,lJ 

Equivalently, (OQ) can be reformulated for high SINR regime and given allocation variables as 
follows: 

L K N / 2 I t \ a n k ,1 
O + in! 



minimize log 2 J] J] J] 

P "' fcJ i= i fe=ln=l \Pri,k,l a n,k,l 

N 

subject to a n> k,iPn,k,i ^ P M». VA;, V/ 



(V) 

N v 7 



n=l 

Note that the numerator in © is a posynomial and the denominator is a monomial, hence © 
is a GP problem in standard form that can be solved optimally through efficient interior point 
methods |fT3l after performing the logarithmic transformation of variables Ifl4ll . However, even 
for small dimensions, it is not recommendable to compute the optimal solution, due to the 
huge computational complexity 0(K LN ) associated with the exhaustive search based subcarrier 
allocation phase. In addition, the GP based power allocation method discussed above has two 
restrictions: high-SINR assumption and centralized time-consuming computations. Due to the 
mentioned facts, there is a need to develop bounds and sub-optimal resource allocation schemes 
for multi-cell OFDMA networks. 
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III. Bounds on the Network Throughput 



A. Lower Bound on the Optimal Network Throughput 

A LB for the optimum multi-cell network throughput can be computed by considering worst 
case ICI. Observing the dependency of ICI on the subcarrier allocation and power allocation 
variables, we assume that each user in each cell is transmitting on each subcarrier with its 
maximum power. A simple LB for the average network throughput C taking the worst case ICI 
into account can be written as follows: 



where £ n>l — J2j=l,j# J2k=l ^fc.raaxflW.fcjJ- 

A tighter LB can be derived by using Algorithm 1 where each subcarrier is allocated to the 
user that maximizes Q n ,k,i where: 



Thus, Q n> k,i is an SINR term for each user k at each subcarrier n in each cell I assuming worst 
case interference. We collect these SINR terms into a vector = [q n ,i,h Qn,2,i----, Qn,K,i] and 
then stack all the vectors in a matrix Q; G IR 7VxX . The resulting allocations based on this criteria 
are then used to compute the LB network throughput using (IB- 
Note that if ^ rit i = 0, than Q n ^,i becomes the marginal rate which is shown to be a near-optimal 
criterion in single cell network scenarios without ICI B2J. Moreover, equal power allocation has 
insignificant performance loss in high SINR regime compared to the optimal water-filling solution 
ED, 01, thus power equalization is implemented in Algorithm 1. For the low SINR regime, we 
can incorporate water-filling rather than equalization in a straightforward manner. 

B. Upper Bound on the Optimal Network Throughput 

Establishing an UB is significantly important in order to calibrate the performance of sub- 
optimal resource allocation schemes with respect to the optimal solution. The UB can be derived 
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Algorithm 1 Computing LB and UB Allocations in Cell / 

1) Input: [Hi], [Ai], [Pi], [G j7 ] where a nA i = 0,p nA i = P k ,m^/N Vfc,Vn 

2) For each user k in cell /, power is divided equally over all of its allocated subcarriers and 
the remaining unallocated subcarriers of the system. 

3) Using [Pi] from step 2, [Hi] and [Gji], compute the matrix Q/ for each cell 

4) Find the (n, k) pair that has the maximum value of Q n ,k,i- Allocate subcarrier n to user k. 

5) Delete the n th subcarrier from the set of unallocated subcarriers. 
If there are still unallocated subcarriers in the system go to step 2. 

else terminate after distributing the maximum power equally at each user over all of its 
assigned subcarriers. 



by ignoring the effect of ICI in all the cells. This can be achieved by substituting £ n j = in 



The allocations computed by Algorithm 1 are near optimal since they are based on a criterion 
which is shown to be near optimal in the context of single cell scenarios GO-H. The average 
network throughput revealed by these allocations could be highly optimistic for multi-cell sce- 
narios. Thus, we can investigate the impact of ICI by simply computing the throughput using 
(OQ) instead of (ITOl ) with these allocations. Computing throughput in this way helps to analyze 
the degradation in the performance when the single cell near-optimal allocations are used in 
multi-cell network scenarios with ICI. 

C. Complexity Analysis 

The (n, k) pair at which the term Q n ^,i becomes maximum is allocated (Step 4), which has 
a complexity of a two dimensional search, i.e., O(KN). However, as soon as a subcarrier is 
assigned, each user updates its power as defined in Algorithm 1. This process iterates until all 
the subcarriers in all the cells are allocated and, thus, the time complexity of Algorithm 1 is 



Algorithm 1, i.e., Q nAt 





(10) 



0(KN 2 ). 
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D. A Motivating Example 

Consider an example with two cells, two users and two subcarriers. Each user can transmit 
with a maximum power of 1 W. Assume Hi= [1 0.9; 0.8 0.7] and H 2 = [1 0.9; 0.8 0.7]. Single cell 
allocation strategies that aim to maximize the local throughput of each cell suggest A 1; P 1 and 
A 2 , P 2 =[1 0; 1]. Computing the UB using CUl results in 1.7655 bps/Hz/cell where a 2 = 1. 
Now, assume the knowledge of interfering link gains at each BS, i.e., Gi 2 = [0.9 0.2; 0.2 0.9] and 
G 2 i= [0.7 0.1; 0.1 0.7]. Computing the throughput again while keeping the single cell allocations 
and taking into account the interfering gains leads to an average network throughput of 1.1137 
bps/Hz/cell. However, better allocations are possible if we consider Ai, Pi and A 2 , P 2 = [0 1; 
1 0] as per the criterion discussed in Section IV which enhances the resulting average network 
throughput to 1.5977 bps/Hz/cell. 

IV. Sub-Optimal Centralized Resource Allocation Schemes 

Considering the high intricacy of implementing the optimal solution, we develop a two-stage 
centralized scheme. In comparison to the centralized scheme presented in IPT21 . the subcarrier 
allocation phase of the developed scheme is iterative and computationally intensive. However, 
the performance is better even without the power allocation phase and, thus, it can provide an 
effective benchmark for low complexity schemes. Compared to 1(121 . we also discuss the power 
allocation phase for the general SINR regime. 

A. Centralized Scheme A 

In the proposed scheme, we split the resource allocation procedure into two phases: subcarrier 
allocation phase and power allocation phase. It is important to note that the subcarrier allocation 
phase involves a power equalization step, thus, it is not totally independent of power allocation. 
Phase I: Subcarrier Allocation 

• Initial Allocation: Firstly, we define the term for the allocation of resources to the users as 
follows: 

Pn,k,lhn,k,l 

Xn,k,l~=Z p (11) 

J j~ 1 J~/~ ! k,maji9n,k,lj 
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This criterion guarantees the selection of the users who possess not only better power- 
gain product but also they offer less interference to the neighbor cells. The denominator 
J2j=i jjki Pk,ma, x g n ,k,ij accounts for the maximum aggregate interference that the k th user 
in cell /, may cause to all cells. Even though this criterion is heuristic, it improves the 
performance compared to the traditional C/I scheme (which gives nearly similar results as 
our lower bound). Once the initial allocations are computed, we can calculate the initial 
throughput of the network C Q using ©. 
• Maximize Throughput Iteratively until Convergence: In this step, we select any cell I and 
subcarrier n arbitrarily and re-perform the allocation at this subcarrier considering the other 
cell allocations fixed, i.e., I n j remains fixed. More explicitly, we compute C n> k,i = log2(l + 
Vn,k,ihn,k,i ^ £ Qr ^ users m i one one anc j se i ect t h e user w hich gives the maximum 

incremental throughput at subcarrier n, i.e., C n — C Q . Note that, in order to compute 
C n ,k.h we need to compute p n ^,i which can be obtained simply by dividing Pfc >max equally 
among all the fixed allocated subcarriers of user k and the new one which is currently under 
observation. 

Once the reallocation is done at subcarrier n, we move to the next subcarrier in cell / 
and so on. As the new allocations are computed for cell I, we calculate the new increased 
network throughput C new and move to another cell j. The whole process is repeated again 
with C Q = C ncw until convergence to a desired accuracy is achieved. 
Phase II: Power Allocation 

Once the subcarrier allocation is done, the optimal powers can then be calculated for the high 
SINR regime or for the general SINR regime through solving a series of GPs using successive 
convex approximation which is a provably convergent heuristic [[T4|. This approach is known to 
compute globally optimal power allocations in many cases. Thus, for given allocations, (OQ) can 
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be formulated for the general SINR regime as follows: 



L K N 




) 




Pn,k,lh n ,k,l + CT 2 + I n ,l 



(12) 




,max; 



Vfc,VZ 



n=l 



Note that the numerator and denominator in (fT2~l) are posynomials and minimizing a ratio between 
two posynomials is referred to be a truly non-convex NP hard intractable problem known as 
complimentary GP. However, this problem can be transformed into GP by letting the denominator 

f(j>) = Pn,k,ihn,k,i + cr 2 + I n ,i = Ya=i Y!k=i u n,k,i(p) and approximating the denominator f(p) 
with a monomial using the arithmetic/geometric mean inequality as follows: 



where s n ^,i = ^yw) • Thus, the problem can be solved by extending the single condensation 
method presented in [TT4ll for multi-cell scenario. The details of centralized scheme A are 
presented in Algorithm 2. 

B. Centralized Scheme A: Complexity Analysis 

The initial allocation phase has a complexity of 0(KN 2 ) which is the same as Algorithm 1. 
Next, we perform a one dimensional search for the user in cell I with maximum incremental 
throughput at subcarrier n. The process is repeated for each subcarrier and cell. Thus, the com- 
putational complexity of this step is O(KNL). Since, the process continues until convergence, 
(i.e., M iterations), the complexity of this step can be written as O(KNLM). Finally, the total 
complexity of subcarrier allocation phase is 0(KN 2 + NKLM). 

The complexity of Phase II is difficult to determine, however, it can be measured in terms of the 
degree of difficulty (DoD) that in turn relies on the number of constraints and variables associated 
with the GP [15J. Since we are dealing with LK power constraints and LKN power variables, the 
total computational complexity of centralized scheme A is 0(KN 2 + NKLM) + DoD(LKN). 
Apparently it seems that implementing centralized GP/successive GP based schemes may not 
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Algorithm 2 Centralized Scheme A 



1) Input: [Hi], [Aj], [PJ, [Gy] where a„ iAi , = 0,p nj/M = P fc , max /iV VA;,Vn 
Subcarrier Allocation (Phase I) 

Initial Allocation: 

2) For each user k in cell /, power is divided equally over all of its allocated subcarriers and 
the remaining unallocated subcarriers of the system. 

3) Using [Pi] from step 2, [Hz] and [Gy], compute Xn,fc,z for every k th user at n th subcarrier 
in cell I. 

4) Find the (n, fc) pair that has the maximum value of Xn,k,i- Allocate subcarrier n to user k. 

5) Delete the n th subcarrier from the set of unallocated subcarriers. 
If there are still unallocated subcarriers in the system go to step 2, 

else terminate after distributing the maximum power at each user over all of its assigned 
subcarriers 

6) Compute C Q 

Maximize Throughput Iteratively until Convergence 

do while(C ncw - C Q > e) 

I = 1, do while I < L, I = I + 1 

n = 1, do while n < N, n = n + 1 

k = 1, do while k < K,k = k + 1 

7) Allocate the subcarrier n to user k. 

8) Compute -p n ,k,i by dividing Pfc,max equally among the allocated subcarriers. 

9) Compute C niM - C Q 
end 

10) Allocate subcarrier n to the user who maximizes C Ut k,i — C Q 
end 

11) Compute C ncw using ©• 

12) C Q Cncw 

end 

Power Allocation (Phase II) 

13) Compute the optimal powers P; in the high SINR regime © given the allocations from 
Phase I. 

14) For general SINR regime, take P^ from step 13 as an initial starting point. 

15) Using Pi, evaluate p n ,k,ih n ,k,i + o" 2 + I n ,i for each allocated user k in cell I at subcarrier n. 

16) Compute the weights s n ^,i as follows: 

_ "n,i,l 

17) Approximate the posynomial using (TOT) . 

18) Solve the approximated GP using any available commercial software iTOl 

19) Go to step 15 using P; of step 18 until convergence. 
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be a good choice for practical implementations. However, in order to reduce the complexity 
and DoD of the power control phase, we have developed the following less complex centralized 
scheme. 

C. Centralized Scheme B 

In this scheme, firstly the subcarriers are allocated in each cell I using the heuristic criterion 
defined in (fTT|) . The allocation of each subcarrier is followed by the power allocation phase 
(based on equalization) as mentioned in the initial allocation phase of Algorithm 2 (i.e., Steps 1 
to 5). Once the subcarrier allocations are finalized, we then compute GP based powers for the 
allocated users at any arbitrarily selected subcarrier n in all cells. Setting the equalization based 
powers p n ,k,i cq as the upper bound on p n ,k,i and considering a high SINR regime, we now define 
the following less complex GP problem with the objective to maximize the throughput at the 
n th subcarrier: 



Clearly, the resulting GP based power of each competing user at subcarrier n in the different 
cells may not succeed in achieving the upper bound, due to the ICI effect. We call this power 
as left-over power. The left-over power can then be distributed equally among the remaining 
allocated subcarriers of the user. The procedure is detailed in Algorithm 3. 

Since at the end of the initial allocation phase, the subcarrier allocations become fixed and 
the total power is distributed equally among the allocated subcarriers of a user, we cannot set 
an upper bound which depicts higher power than the previously allocated power. If we do so, 
this may cause power reduction or even no power at some other allocated subcarrier of that user 
in order to maintain the total power constraint. Thus, this may results in an invalid subcarrier 
allocation. 

Next follows an example which demonstrates the significance of GP as well as centralized 
scheme B over equal power allocation. Consider H 1 ,H 2 = [0.30 0.25; 0.04 0.15] x 10~ 9 , 



minimize log 2 

Pn,k,l 




1 = 1 



(14) 



subject to p nAl < p n ,k,i cq , V/ 
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Algorithm 3 Centralized Scheme B 

1) Repeat Steps 1 to 5 of Algorithm 2, i.e., initial allocation phase. 
n = 1, do while n < N, n = n + 1 

2) Compute the GP based powers p Hi k,i of the allocated users at any subcarrier n considering 
a high SINR regime using (fl4l) . 

3) For each user allocated in a cell I at any subcarrier n, divide the left-over power equally 
among the remaining allocated subcarriers of the user. 

4) Remove the subcarrier n from the set of unallocated subcarriers. end 



da = [0.06 0.05; 0.16 0.06] x KT 11 and G 2 i=[0.14 0.69; 0.76 0.1935] x 1CT 11 . The equal 
power allocations dictate Pi = [0 0.5; 0.5] and P 2 = [0.5 0; 0.5 0] which leads to an 
average network throughput of 11.8392 bps/Hz/cell. However, computing the GP based powers 
results in P 1 = [0 0.53; 0.47] and P 2 = [0.38 0; 0.62 0] which lead to a maximum 
average network throughput of 17.2734 bps/Hz/cell. 

D. Centralized Scheme B: Complexity Analysis 

The initial allocation phase has a complexity of 0(KN 2 ) which is the same as Algorithm 1. 
Since (fT4l) has L constraints and variables, the complexity of the power control phase is signifi- 
cantly reduced. Although this procedure restricts the degree of freedom offered by GP, numerical 
results show that the network throughput remains comparable with reduced complexity. 

V. Distributed Resource Allocation Scheme 

In the centralized strategy, we assume that Xn,k,i is known, i.e., every BS knows the interfering 
gains offered by its users to the neighboring BSs. The interfering gains are based on path loss, 
shadowing and fading. Assuming the knowledge of local user positions at each BS, the path 
loss of local users toward the first tier of interfering cells can be determined, however, the 
knowledge of shadowing and fading gains is difficult to assume in practical scenarios. Thus, in 
the distributed approach, we compute our results without using the knowledge of shadowing and 
fading interfering gains. 

Each BS performs the subcarrier allocations without taking ICI into account. In other words 
we compute single cell near optimal allocations using Algorithm 1 . The allocation decisions are 
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locally made at each BS and do not need collaboration. Once the allocations are decided, each 
cell shares them with all other interfering cells. The GP based optimal powers in © can then 
be evaluated in a distributed way using dual decomposition methods by first performing the log 
transformation of the variables, i.e., \np n k j = p n ^,i and lnp n)fej - = p n ,kj, then adding auxiliary 
variable \nz n ^j = z n ^j where z n ^ = p n ,k,j m order to transfer the coupling in the objective to 
coupling in the constraints lfT4ll . For given allocations, the problem in © can thus be written in 
a distributed way as follows: 



L N 



minimize log 2 

l=l n=l \ ' ' 

N 

subject to ePn ' k ' < p fc,max, Vfc, VZ (l5) 

n=l 

Z n ,lj = Pn,k,j, Vn, V/ 

Since the computational complexity of (fT5l) is high as it has LK power constraints and LKN 
variables, we present the dual decomposition of (fl4l) which is more suitable for practical scenarios 
and has a lower computational complexity. Moreover, the objective function in <fT5l) not only 
depends on the powers of local users p n ,k,i but also on the power of users sharing the same 
subcarrier in neighboring cells p n ,k,j- Thus, in order to minimize the objective in (fT5l) . each 
BS requires the knowledge of interfering gains and interfering transmit powers, that may lead 
to significant overhead to exchange control information. Thus, in order to obtain a practical 
distributed solution, we keep a local copy of each of the effective received powers i.e., z n ij = 
9n,k,jiPn,k,j D31- CH1> can then be formulated in a distributed way as follows: 

minimize > log 2 ^ — — 

Zn,lj,Pn,k,l ^ y e Pn ' k ' l K,k,l 

subject to e pn ' k ' 1 < P n , fc ,z,e q , V/ 

%n,lj 9n,k,jl Pn,k,j 



° 2 + Y!j=l,j^l 9n,k,jl e ' 



n,lj 
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The partial lagrange C{$> n ,k,h Zn,ij, ^> Vn,ij) f° r (f!6l ) can then be written explicitly as follows: 

Yl l0g 2 e^Zh^ 1 + S Vnjj {Zn,lj ~ gn,kjl ~ Pn,k,j) + J2 A ' ~ P «,M,eq) 

1=1 \ ' n,fc,/ / 1=1 j=l,i+l 1=1 

(17) 

Eq. (TT71 ) can be decomposed into L sub-problems with local variables p nt k,h z ni ij, M an d coupling 
variable r] ni ij. The simple lagrangian Ci for each cell I can then be written as follows: 



( a 2 + ]T - e ,u \ / \ 

A = l0g 2 + Vn,ljZn,lj ~ ^ ^ n ' fc >' + A ' ^ n ' fc '' ~ P '«,M,eq) 

V 6 ^ / i=l,j+l \j=l,3+l J 

(18) 

where Xi is the lagrange multiplier for the inequality constraints and r\ n ^ 3 - are the consistency 
prices. Thus, the minimization of (fl"8l) with respect to the local variables can be done in a 



distributed way at all BSs. At every iteration, each cell / receives the term [ ^2 r] n ji J by 
message passing and minimizes the local Lagrangian (1181) with respect to the local variables 
Pn,k,h z n ,ij, h subject to the local constraints. In order to obtain r\ n ^ the following master lagrange 
dual problem has to be solved: 

L 



maximize minimize Ci (19) 
A simple way to solve (fl9l) is to use the following subgradient update for the consistency prices: 

Vn,lj{t + 1) = Vn,lj{t) + (S/t) (Znjj - log 2 Pn,k,j9n,k,jl) (20) 



In summary, each BS minimizes (1181) in parallel with respect to the local variables after receiving 
the term Ylj=i j+iVn,ji- Each BS then estimates the received interference z n> ij from each cell 
and update the local consistency prices using (1201) . Finally, each BS broadcast them by message 
passing to all BSs. Note that 5 in (1201) represents the step size and is non-negative. 



VI. Performance Evaluation 

A cellular OFDMA network is considered where the radius of each cell is assumed to be 
1 km. The maximum user transmit power is considered to be 1 W. The channel gain is defined 
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as follows: 

h n>k ,i = (-122 - 10 7 log 10 4,«) - Af(0, a 2 ) + 10 \og w F nAl d > d iei 
h n>k>l = (-122 - IO7 log 10 d rcf ) - JV(0, a 2 ) + 10 \og w F nAl d < d Te{ 

where d TC { is the reference distance and is set equal to 0.05 km, dk,i is the distance of the k th 
user from the I th BS. The first term denotes the path loss where 7 is the path loss exponent 
and is set equal to 3. The second term represents log-normal shadowing with a mean of dB 
and a standard deviation of 8 dB. The last factor F n ^,i corresponds to Rayleigh fading. The 
bandwidth of the system is assumed to be 20 MHz with a noise power spectral density of 
8.6455 x 10~ 15 W/Hz at each receiver. The channel conditions are assumed to be fixed during 
a frame. The interfering gains from the j th interfering cell to the cell of interest I are computed 
as follows: 

9n,kji = (-122 - 10 7 logio4,0 - -A/"(0. ° 2 ) + 10 logio^n,fc,i 

where d^i is the distance between the k th user in the interfering cell j and the I th BS. We 
consider the following two simulation scenarios: 

• Scenario A: Users are equidistant from the BS and placed at equally spaced angles. 

• Scenario B: Users are assumed to be uniformly distributed across the whole cellular area. 
In Table 1, we compare the performance and complexity of the centralized and distributed 

schemes with the derived bounds and the optimal solution in high SINR regime. The optimal 
solution is computed by the exhaustive search based subcarrier allocation phase detailed in 
Section EE. All users are placed at equal distance d from the BS and at equally spaced angles 
(i.e., scenario A). The results are taken after averaging over 100 channel realizations. The 
simulation results show that the performance gap between the benchmark centralized scheme 
A (with power control) and the optimal solution is negligible compared to the low complexity 
centralized and distributed schemes. However, this observation may not remain valid for bigger 
network scenarios. Moreover, as d increases the degradation of the average network throughput 
is evident. 
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In Fig. 1, Fig. 2 and Fig. 3, we present the performance of the centralized scheme A, centralized 
scheme B and the distributed scheme for two cells, four cells and seven cells, respectively. The 
results have been taken after averaging over 10,000 channel realizations and are shown for 
both simulation scenarios. The performance of all schemes is calibrated using the established 
upper and lower bounds. Since the centralized scheme A has computationally exhaustive power 
allocation phase, the results are presented for the subcarrier allocation phase of Algorithm 2 
only. However, it can be observed that the scheme still has the capability to serve as a suitable 
benchmark for the developed low complexity schemes. In order to highlight the significance of 
the less complex GP problem defined in (fl4"l) . we also present the performance of the centralized 
scheme B without power control. 

For the two cell scenario, the performance gap between the centralized schemes is negligible 
and they give nearly similar results. However, as the number of cells increases the performance 
gain of the centralized scheme A is evident over all schemes even without power control. 
Moreover, it is also important to note the significant degradation in the performance of centralized 
scheme B without power control phase. This degradation is found to be increasing with the 
increase in number of cells. It is also worth to mention here that the proposed less complex GP 
problem (fl4l) can be implemented in a distributed way using the techniques explained in lfT4ll 
and can be used with any set of subcarrier allocations. Thus, in the distributed approach we 
use the near optimal single cell allocations in conjunction with the less complex GP problem 
(fl4"l) . The significance of the power control phase can be observed easily from the results which 
becomes more evident for high number of cells . 

Moreover, the presented results depict the reduction in the average network throughput as the 
number of interfering cells increases. The performance gap of the proposed schemes increases 
with respect to the evaluated UB. Even though the UB is not tight and reflects an over optimistic 
average network throughput, it provides an idea on the performance gap between the proposed 
schemes and the optimal solution. 

In Fig. 4, we assume that the users in each cell are placed at equally spaced angles from 
to 27r. The performance evaluation of all schemes has been done by changing the user positions 
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from cell center to the edge of the cell. It is observed that the performance gap increases between 
the centralized and distributed schemes as users approach the boundaries of the cell. 

VII. Conclusion 

In this paper, we developed an upper bound and a lower bound to the optimal average network 
throughput in multi-cell uplink OFDMA networks. We also investigated the severe effects of ICI 
on the performance of a single cell near optimal resource allocation scheme. Moreover, we 
proposed a benchmark centralized scheme which is useful to study the performance gap of 
the low complexity centralized and distributed resource allocation schemes developed later with 
respect to the optimal solution. All schemes are compared to the exhaustive search based optimal 
solution and derived upper and lower bounds for various scenarios. The schemes are evaluated 
and compared in terms of network throughput and computational complexity. 
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TABLE I 

Average network throughput (in bps/Hz/cell) of the derived bounds, centralized and distributed 
schemes for l=2 cells and n = 6 subcarriers/cell 





K=2 Users 


K=4 Users 


A'=6 Users 


Computation Complexity 




d=0.5km 


<i=0.9km 


d=0.5km 


<i=0.9km 


d=0.5km 


<i=0.9km 






UB 


44.2642 


33.1294 


55.7414 


42.8390 


60.6901 


49.6214 


0(KN 2 ) 




Optimal 


37.1168 


29.8642 


47.9975 


35.6520 


52.1299 


41.0121 


0(K NL ) - 


h DoD(LKN) 


Centralized A 


36.8061 


28.6973 


46.4765 


34.0713 


51.2868 


40.5845 


0(KN 2 4 


- NKLM) + DoD(LKN) 


Centralized B 


36.4755 


27.0352 


45.6239 


33.4280 


49.7971 


38.7237 


0(KN 2 ) 


^DoD(L) 


Distributed 


35.3623 


25.9976 


43.5918 


31.9231 


48.8887 


38.0050 


0{KN 2 ) 


^DoD(L) 


LB 


35.0966 


25.8635 


42.5509 


31.0261 


48.1571 


37.7996 


0{KN 2 ) 
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Fig. 1. Comparison of all proposed schemes for L=2 cells, (a) Scenario A (b) Scenario B: Users are placed at 0.9 km from BS 
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Fig. 2. Comparison of all proposed schemes for L=4 cells, (a) Scenario A (b) Scenario B: Users are placed at 0.9 km from BS 
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Fig. 3. Comparison of all proposed schemes for L=l cells, (a) Scenario A (b) Scenario B: Users are placed at 0.9 km from BS 
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Fig. 4. Comparison of all proposed schemes for L=l cells, iv = 100 users from cell center to cell edge 
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